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ABSTRACT 


a : F : : 
Traverse is a method of surveying in which a sequence of 


lengths and directions of lines Letween points on the Earth 
are measured and used in determining positions of the 
points. This method is one of several used to find the 
accurate geodetic positions which various agencies use. 
Traversing is a convenient, rapid method for establishing 
horizontal control. 

The theoretical background is provided here to explain 
the method of traverse station position computations and 
adjustments in the Universal Transverse Mercator grid coor- 
dinates. Closed traverse station positions were computei 
and adjusted using J the APPLOTi Eat Method and by the Least 
Syuares Methodf *the ‘adjusted coordinates of both methods 
were transformed from the Universal Tcansverse Mercator grid 
coordinates to geographic coordinates and compared with the 
coordinates ynick were adjusted by the U.S. National Ocean 
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4 I. INTRODUCTION , 


Hydrographic surveying includes many branches of science 
for the purpose of the production of nautical charts 


: 
specially designed for use by the mariner. The determina- a 
} tion of position in hydrographic surveying is as important 

J as the measurament of depth. Before determining an accurate 

hydrographic position, accurate geodetic positions! for 

l shore control must be established. There are many methods » 
available to establish geodetic control in the survey area. 

These methods are: 

z 1. Trianyulation 

; 2. MTrilateration » 
| 3. Traverse 

; 4. Intersection 

: 5. Resection 

4 The process of making a proper nautical chart consists, > 
: first of all, in setting up a framework of marks on the 

ground. Before 1950 the main framework of a geodetic survey 
almost always consisted of triangulation, which was replaced 


by traverse if the topography made triangulation impracti- 
cable. During the last decade, the introduction of elec- 
tronic distance measuriny (EDM) equipaent has made both 
trilateration and traverse economical, and an acceptable 
substitute for triangulation. In fact, it appears probable ) 
that these new methods will replace triangulation as the 

fain framework for new geodetic surveys [Ref. 1]. It is 

evident that triangulation and traversé are the main methods 

used for establishing control. There seems to be no - 


1A position of a point on the surface of the Earth. 
expressct in terms of geod tie ent euie. geodetic longitude oe 
ard geodetic height. geodetic position implies an adopter oe] 
geodetic datun. ‘ 
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agreement among the various agencies as to which of these 

two methods is mostly used. The U.S. National OScear Service 

(JOS) does the majority of its horizontal control surveys ; 
for hydrography with traverse (about 90%) [RKef. 2]. The 

main factor for the selection of one or the other method 

depends on thea geographical configuration of the survey on 
project area and the availability of good EDM equipment. 

Traversing is a convenient, rapid method for establishing 

horizontal control. It is particularly useful in densely 

built areas, when the coastline tends to be even, along a 

railroad track, and in heavily forested areas where lengths 

of sight are short so that neither triangulation nor 

trilateration is suitable. 

The objectives of this thesis are to show methods of 
traverse station position computation and a comparison of 
methods of least squares adjustment. Closed traverse 
station positions were computed and adjusted using the 
Approximate Method and by the Least Squares Method in the 
Universal Transverse Mercator (UTM) grid coordinates. The ‘ 
computer programs were written to calculate the adjusted 
traverse station positions. Test data included those data oe 
obtained during the Geodetic Survey Field Experience course bes 
at the Naval Postgraduate School (NPS) in October 1983. 

The results of comparative computations are shown to 
more Significant figures, in this thesis, than are normally 
considered desirable in production work. The same observed 
data are used with each of three computational methods. It 
is important to recognize that what is being compared is not 
observational precision but computational precision. MYence, 
it is considered necessary for a rigorous comparison of : 
computational precision, including round off error, to show : 


resuits to several more decimal places than is justifiable 


based on okservational precision alone. 


eet ee 


II. TRAVERSE = 


A. GENERAL 


Traverse is a method of surveying in which a series of 
straight lines connect successive established points along 
the route of a survey. An angular measurement is taken 
using a theodolite at each point where the traverse changes 
direction. Distances along the line between successive 
traverse points are determined by EDM equipment. The points 
defining the ends of the traverse lines are called turning 
points, traverse points, or traverse stations. Each 
straight section of a traverse is called a ley or a traverse 
line. 

A closed traverse originates at a point of known posi- 
tion and is closed or another point of krown horizontal 
position. Traverse I-1-2-3-F originates at point I with a 
backsight along line IA of known azimuth and closes on point Cae 
F, with a foresight along line FB also of known azimuth 
(Tigure 2.1). This type of traverse is preferable to ali rea 


Figure 2.1 Closed Traverse. 


others since computational checks are possible which allow 
detection of systematic errors? in both distance and 
direction. 

A closed-loop {closed polygon) traverse is a special 
case of a closed traverse in which the originating and 
terminating points are the same point with a known hori- 
zontal position. Traverse I-1-2-3-4-I, originates and 
terminates on point I (Figure 2.2). This type of traverse 
permits an internal check on the angles, but there is no 
check on the linear measurements. Therefore, there is a 
possibility that an error proportional to distance may occur 
and not ke detected. 


Figure 2.2 Closed-loop Traverse. 


B. ANGLE AND DIRECTION MEASUREMENT 


Angles and directions may he defined by means of bear- 
ings, azimuths, deflection angles, angles to the right, or 
interior angles. These quantities are said to be observed 


2Systematic errors may be caused By fous instrumerts 
oc factors such as temperature or humijity changes which 
affect the pecformance of measuring instruments. 
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when obtained directiy in the field and calculated when 
obtained indirectly by computation. 

A theodolite is an instrument desijnec to observe hori- , 
zontal directions and measure vertical angles. It consists 
of a telescope mounted to rotate vertically on a horizontal 


axis supportei by a pair of vertical standards attached toa on 
revolvable circular plate containing a graduated circle for > 
observing horizontal directions. Another graduated arc is 

attached to one standard so that vertical angles can be 

measured. 

Repeating and direction theodolites have features that ’ 
are common to poth types of instruments. Repeating theodol- 
ites are read directly to 20" or 01! and by estimation to 
one-tenth the corresponding direct reading. Direction theo- 
dolites are usually read directly to 01" and can be esti- —° 
mated to tenths of seconds [Ref. 3, p. 215]. In general, 
direction theodolites are more precise than are repeating 
theodolites. 

The direction theodolite observes directions only and » 
angles are computed by subtracting one direction from 
another. Assume that the horizontal angle AIB (Figure 2. 3) 
is to be measured with a direction theodolite. The theo- 


dolite is set over point I, leveled and centered, anda 


sight is taken on poirt A. The Lorizontal circle is then 
viewed through the optical-viewing system and the circle 
reading is observed and recorded. Assume the reading is ce 
459 02" 4O". The telescope is then sighted on point B. The _ 
hocizontal circle is then viewed through the optical-viewing 
system and the circle reading is observed and recordei. 
Assume the reading is 1249 11' 59", These two observations 
constitute directions which have a common reference direc- . e 
tion that is completely arbitrary. The clockwise horizontal 
angle is 

i = (1260 171% 59") - (45° 02* 40") = 790 09% 19" 


A 


/ 


Figure 2.3 Horizontal Angle. 


C. DISTANCE MEASUREMENTS 


There are several methods of determining listance, the 
choice of which depends on the accuracy required and the 
cost. For example, tacheometry, taping, and ED“ equipment 
can be used. The general availability of EDM eguipment has 
practically eliminated the use of taping for the measurement 
of traverse lengths. Accuracies are comparable, or 
superior, to those obtained with invar tapes. 

Distances measured using EDM equipment are subject to 
errors arising from the instrumental components, calibration 
of the equipment, inaccuracies in the meteorological data, 
elevation discrepancies, and centering of the instruments or 
reflectors. The reduction of measured distances involves 
converting the slope distance to a horizontal distance, 
converting the chord distance to an equivalent arc distance, 
and reducing the arc distance to the ellipsoid. The reduc- 
tion of slope distance to horizontal distance is necessary 
to compensate for the difference in elevation of the end 


points of the measured line. The horizontal distance is 


15 


reduced to the geodetic distarce Ly applying a sea level 


corrector ard a chord-arc corrector to the horizontal 
distance (Ref. 4, pp. 124-125]. 


D. ACCURACY 


In general, the accuracy of a traverse is judged on the 
basis of the resultant error of closure of the traverse. 
This resultant closing error is a function of the accuracies 
in the measurement of directions and distances. The classi- 
fication, standards of accuracy, and general specifications 
for horizontal control have been prepared by the Federal 
Geodetic Control Committee (FGCC) { Ref. 5] and have been 
reviewed by the American Society of Civil Engineers, The 
American Congress on Surveying and Mapping, and the American 
Geophysical Union (Table I). The third-order traverse class 
I and class II are of particular interest because these are 
the orders of accuracy the hydrographer is usually reguired 
to accomplish. 


E. DATA ACQUISITION 


The data for this thesis were acquired at Moss Landing, 
California, by NPS Hydrographic students during the Geodetic 
Survey Field Experience course in Ictober 1983. All of the 
known station positions and azimuths (Table II) were 
adjusted by the Coast and Geodetic Survey by third-order 
methods [Ref. 6]. 


The distances were observed by using a Kern DM102 anda 
Tellurometer MRAS. The Kern DS102 is an electro-optical 
distance meter. The measuring accuracy was #(5 mm + 5 ppm). 
The Tellurometer MRA5 is a microwave instrument. Precision 
in terms of probable error, in the temperature range of 
-329C to +440C, is +(5 cm + 100 ppm) for a Single 
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TABLE II 
Coordinates of Known Stations and Azinuths 


= oe Naan lag Sa i ea 
. _, UTM grad coordinates . | 
Station as La ing Grid Easting | 
Mm. 7 me 
Moss 2 ; 4,072,555. 85206 608,279.044n4 
Holn 4,079,258.31754 612, 238.85256 
Grid azimuth clockwise ' 

| from North Azimuth | 
(0) e " | 
From Moss 2 to Pipher stations 100 16 23.778 | 
From Foln to Moran stations 136 33 26.334 | 


—— a ee ee = Se 


determination. The distances were observed in the field, 
cocrected by temperature and pressure for propagation error. 
Unfortunately, meteorological data are generally acquired 
only at the end points of a measured line. Using the mean 
of these meteorological values only approximates the actual 
conditions of the entire measured line ani does not 
completely correct for errors in the propagation velocity of 
electromagnetic radiation. By applying the elevation, sea 
ievel, ard scale factor corrections they were reduced to UTM 
grid distances [Ref. 4, pp. 124-125]. The OTM grid 
distances and standard deviations of the distances were 
determined in meters (Table III). 


The angles were observed by using a Wild T-2 theodolite. 
To ensure the correctness of the beginning and ending 
azimuths, a check azimuth to a second station of known posi- 
tion was observed. The angles were observed at stations 
Moss 2, Mossback, Dune Temp, and Holm (Table Iv). All 
observations were made by NPS students and conform to 
specification for a third-order class I traverse. 
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TABLE II1I 


Grid Distances and Standard Deviations 


"4 
| 


P Standard 
Grid distances between Distances deviatior 
stations ( mo. ) ( m. 
Moss 2 and Mossback 1,424.004 0.001 | 
Mosskack and Dune Temp 365.744 0.001 | 
Dune Temp and Holn 6,476.271 0.003 | 
—_ pleases ae SERS a A SS eR ETS 
TABLE IV 
The Observed Angles and Standard Deviations 
ee eee ee eee Te TE ae ee eee ee ee ee ee ee 
| 
Backsight Center Foresight Observed Standard | 
station station station angles deviation | 
, 0 t CL] " 
Pipher Moss 2 Mossback 246 05 43.200 01.984 
Moss 2 Mossback Dune Temp 222 51 98.600 01.405 
Mossback Dune Temp Holm 190 15 902.600 91.203 


Dune Temp 


Holm 


| 
Moran 277 O5S 17.000 01.614 | 


ae 


aos. 


III. TRAVERSE COMPUTATIONS AND ADJUSTMENTS 


A. INITIAL TRAVERSE COMPUTATIONS 


A traverse (Figure 3.1) originates at statior 1 (known 
position) and terminates at station 4 (also known position) 
Table V). To. compute the forward azimutk (Table VI) of an 
unknown leg, the angle is added to the back azimuth of the 


previous leg (Equation 3.1). 


! 
az; = Az. + 2 oj - (i-1) 180° (3.1) 
Where: i = the number of legs, 
Az, = the forward azimuth of an unknown leg, and 
Az = the fixed initial azimuth. 


For the computation of UTM grid coordinates, let OE; 
and Q\N, be designated as the departure and latitude for leg 
i (Figure 3.1) so that the general formulas to compute the 
departure and latitude are 


OE; = a; sin Az; (3. 2) 
ON, = ad; cos Az, (3. 3) 
where Se Te ee 2g 3 gk 1S ei iDy 

n = the number of the observed argles. 


The algebraic signs of the departure and latitude for a 
traverse leg depend on the siyns of the sine and cosine of 
the azimuth of that leg. The algebraic sign of the depar- 
tuce and latitude is determined by the following rules: 

1. UTM grid azimuths are refered to grid north. 
2. For azimuths between 0° and 180°, the departure is 
plus; for all other azimuths, the departure is minus. 


3. For azimuths between 90° and 2709, the latitude is 


minus; for all other azimuths, the latitude is plus. 


etn Me 


[eee eee ee eer ee en aay = 


Figure 3.1 A Closed Traverse. 


The calculators and computers with internal routines for 
trigonometric functions yield the proper sign automatically 
given the azimuth. The major portion of traverse computa- 
tions consists of calculating departures and latitudes for 
successive legs and cumulating the values to determine the 


coordinates for consecutive traverse stations. 
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TABLE V 
Data for Initial Traverse Computations 


ee ee ee ee ee ne ee ere ne Ce ee a ee eee Ee ee , 
Leg re Angle | ; 
ne 2 
0 t " z 
t= 2 ad, = 1,424,004 A, = 246 65 43.200 
2.=.3 d, = 365.744 QO, = 222 51 08.600 | 
3-4 d, = 6,476.271 QO, = 190 15 02.600 | 
The initial UTM grid azimuth Az_ = 100 16 23.778 | 
{ 
: UTM grid Northing N, = 4,072,555.85206 mn. } 
Station 1 : . | 
OTM grid Easting E, = 608,279.048404 on. 
a Sian a | 
The coordinates of station j {j = i+ 1) are E; and 
Nie the coordinates of station j are 
BE, = B + OF, (3.4) F. 
No = 8, + ON, (3.5) 


J { 


where i= 1,2,3,... NR. 


The values of departures, latitudes, and coordinates are 
Shown in Table VII and were computed from the data in Tables 
Vand VI by application of Equations (3.2) through (3.5). 


B. COMPUTATION OF DISCREPANCIES 


A closed traverse at Moss Landing (Figure 3.1) origi- 
nated at station 1 (Moss 2) ard closed at station 4 (Holm). 
The closing errors for a traverse are caused by observation 
errors in the observed angles and the measured distances. 
The closing errors may be computed by applying Equations 
(3.1) through (3.6). 
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TABLE VI 
The Azimuth Calculation 


| aes SS ee ee eR ge 1 , 
Forward Back Observed | 
Azimuth Azimuth Angle ; ; a 
AZe 100 16 23.778 : 
a, 246 05 43.290 p 
| Az, 346 22 06.978 | 
- 189 00 00.000 | 
| AZ, _, 166 22 96.978 | 
a, 222. 51 —08. 600 , 
AZo 389 13 15.578 | 
- 180 00 00.000 
| AZ3_5 209 13 15.578 | 7 
a, 190 15 02.600 p 
Az, 399 28 18.178 | 
- 180 00 00.000 | 
AZ, _, 219 28 18.178 | -- 
Sa EL ee am ee ee ED ae ae ee : 
_ The angular error of closure may be computed by applying . 
Equation (3.1) becomes , 
n 
Az. + 2a,- (n-1) 1800 - Az, = fy (3-6) 
’ 
where Az, is the fixed closing azimuth and W, is the 
angular error of closure. AZ is the grid azimuth fron 
stations Holm to Moran (AZ, = 1369 33" 26.334") (Takle ITI). 
The computed closing azimuth is egual to AZ, _3 plus the . ’ 
observed angle at station Holm = (219° 28 18.178") + 
(277° O5* 17") = 4969 33" 35.178" = 1369 33% 35.178". From 


Eyguation (3.6), the angular error of closure is es 
(136° 33* 35.178") =~ (1369 33" 26.334") = + 08.844", ' 
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TABLE VII 


Initial Traverse Computations 


ee en ee ee eg ee en wad te We tg Pom 
Values in meters Values in meters | 

N, 4,072,555.85206 E, 608,279.04404 | ora 
ON, 1,383.89267 QE, - 335. 60166 | 
No 4,073,939.74473 Es 607,943.44238 
ON, 319.20061 E5 178.54871 | 
N5 4,074, 258.94534 Es 608,121.99109 
ON, 4,999.282384 QE, 4,116.94755 | 
Ny 4,079, 258.22818 E, 612,238.93864 | 
Ue erate ieee ed 


To compute the errors in closure in position for a 


closed traverse. Equations (3.2) to (3.5) are applied to 


obtain: _ 
n-1 
W, = 2 Oe - (EB, -£,) (3.7) 
n-1 We ee 
Ww, = 24K -(N,- N,) (3. 8) wee 
or os sae 
We BY BRS Be (3.7a) 
W, = N, - NL (3.8a) 
Where: n = the number of observed angles, 
Wo , Wz = precalculated discrepancies, 
E, 5 N, = the fixed initial coordinates, 
E, , Ny = the fixed closing coordinates, and 
E, « N, = the computed closing coordinates. 
The result of the errors in closure in position of a oR 


traverse is 
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Qa = ( W2 4 W212 (3.9) 


The computed closing coordinates of this traverse are 


E, = 612,238.93864 m and N, = 4,079,258.22818 a. The fixed 


closing coordinates are E, = 612,238.85256 m and NL = 


4,079,258.31754 m. The precalculated discrepancies can be e 
computed by using Equations (3.7a) and (3.8a). The results 
are W, = 0.08608 m and W, = -0.08936 m. By applying 

Eguation (3.9), Qd is 0.12408 m. From Table III, the total 

measured distance is 8,266.019 m. The ratio of the distance e 
error of closure to the total distance is 0.12408/8,266.019 

or 1 part in 66,618. The ratio is an indication of the 

goodness of this traverse. 


of closure of this traverse is +8.844", which corresponds to 


o 

C. ADJUSTMENT OF A TRAVERSE BY AN APPROXIMATE METHOD ae 
The angular error of closure of traverse may be distrib- oe 
ute1 equally among the observed anyles. The angular error = 


a correction for each angle of -2.211". The observed angles 
were correctei by this value. The adjusted azimuth was then 
computed. The departures, latitudes, and coordinates were 


recomputed Fy using the adjusted azimuth (Table VIII). The 
closure corrections in E and N coordinates are +0.09635 m 
and -0.04326 m, respectively. The resultant closure is 
0.10562 m. The ratio of the distance error of closure to 


the total distance is 0.10562/8,266.019 or 1 part in 78,262. 2 
The adjustments of a closed traverse by the Approximate 

Method are completed by the compass rule which proportions one 

the errors in E and N coordinates according to the distance a 


of the course [Ref. 4, pe. 354]. The corrections are applied 
to the departures and latitudes prior to computation of 
coordinates. 
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TABLE VIII 


Initial Traverse Cosputet sous. (Approximate Method with 
Adjusted Azimuth) 


ea ee 
Observed angle Correction Adjusted azimuth 
0 e " 
AZ. 100 16 23.778 
GQ, 2469 05" 43.2" - 2.211" 246 05 40.989 
AZ, 346 22 94.767 
- 180 00 00.000 
AZ>., 166 22 04.767 
a, 2229-54") 0826" «= 22! 222 51 06.3895 
AZo 389 13 11.156 


- 180 00 90.000 
209 13 11.156 
a, 190° 15* 02.6" - 2.211" 190 15 00.389 
399 28 11.545 
- 180 00 00.000 


a ee re a ce Sr 
~ 
N 
wW 


Az 219 28 11.545 
QO, 277° 05" 17.0" = 2.2118 277 (05 14.789 
AZgo 136 33 26-334 
Values in meters Values in meters 
N, 4,072, 555. 85206 z, 608,279. 04404 
ON, 1,383.88907 AE, -_335.61649 
; No 4,073,939.74113 Eo 607,943. 42755 
ON, 319.20444 QE, 178.54187 
N, 4,074, 258.94557 E, 608,121.96942 
ON, 4,999.41523 QE, 4,116. 78679 
4,079, 258. 36080 E, 612,238.75621 
7 42079, 258. 31754 Er  612,238.85256 
an = 0.04326 dE - 0.09635 


4 
Ww 
Ne es a ee ae a eS ne A cc Ea SS yA wanes SO a ny ee Eo SS aa te A ty sete a yas 


eee ae ae re ee ee ree te ee ee a a a ee ae ee a 


7 SE. = (dE/D). 4; (3.10) 
: 
Oy; = ( aN 7D) ~ a, (3.11) 
Where: SE, = correction to QAE;, 
5N. = correction to ON,, 
dE = total closure correction in the © coordinate, 
agN = total closure correction in the N coordinate, 
and 
D = total distance. 


The adjusted E and N coordinates (Iable IX) were 
different from the fixed E and N coordinates (Table II) by 
0.00001 m due to round off error. The calculations and 
adjustments were illustrated in sections A, B, and C by 
using the hand calculator (TI-59) and rounded off at 5 
decimal places. The computer program was written to calcu- 
late the adjusted traverse station positions by the 
Approximate Method by using 16 decimal places (Appendix A). 
The approximate traverse adjustment is based on an assumed 
condition that the angular precision equals the precision in 
linear distance. 


D. ADJUSTMENT OF A TRAVERSE BY LEAST SQUARES METHODS 


The method of least squares provides a rigorous 
adjustment and best estimates for positions of all traverse 
Stations. The Least Squares Method is used to 
Simultaneously eliminate closing errors in azimuths and 
coordinates of traverses. 


1. Zhe Principle of Least Squares 


The fundamental condition of the least squares tech- 
nigue in surveying requires that the sum of the squares of 
the residuals be minimized. A residual is defined as the 
difference between the true and observed values. In making 
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TABLE IX 
Adjusted Coordinates (Approximate Method) 


<< SO eee 


cs 


| 
| 
Departure Correction Grid Values in meters 
E, 608,279.04404 
QE, ~335.61649 + 0.01659 - 335.59990 | 
E> 607,943. 44414 | 
QE, 178.54187 + 0.00426 178.5 4613 | 
E. 608,121.99027 | 
QE, 4,116.78679 + 0.97549 4,116. 86228 | 
E, 612, 238.85255 | 
| 
Latitude Correction Grid Values in meters | 
| N, 4,972,555.85206 | 
QN, 1,383.88907 - 0.00745 1, 383. 88162 J 
N, 4,073,939.73368 | 
ON, 319.20444 - 0.00191 319. 20253 
Ny 4,074,258.93621 
ON, 4,999.41523 - 0.03389 4,999. 38134 | 
N, 4,079, 258.31755 | 
a eee 


physical measurements, the true values can never be deter- 
Mined. The least squares principle establishes a criterion 
for obtaining the best estimates of the true values. 

If the best estimates of the true values are stated 
by x, 


; and observed values by Xie the residuals are expressed 


as 
v; = x; = x; 
The fundamental condition of least squares, for 


uncorrelated observations with equal precision are expressed 


as 


n 

it ¥;,)? = (v,)? * (v,)? ¢ (¥,)2 + 2 « - + (¥,)? = Mininua 
i:1 
or in matrix forn 


viy = Minimun (3. 12) 


In general, the observed values are of unequal 
precision. The observed value of high precision has a small 
vaciance. Conversely, a low precision of the observed value 
has a large variance. Since the value of the variance goes 
in opposite directior to that of the precision, the observa- 
tion is assigned a value called weight corresponding toa 
quality that is inversely proportional to the observation's 
variance. 

For uncorrelated measurements x,, Xoo Xan ee 


2 2 
Xn . With variances 0,. Ox o2 erg? os JO: respectively, the 


3° n? 
weights of these uncorrelated measurements are, 


2 2 2 2 2 2 
P,, = 0.7 6, Poa 0,/ O56 Po = 07 om «2 @ 
2 2 
Peas C7 6, (3. 13) 


2 
where 6, is the proportionally constant of an oftservation 
of unit weight [Ref. 7, p. 67]. These weights may be 
collected into a corresponding diagonal P matrix, called the 


weight matrix: 


pP,9 0. - 0 
¢) 52? 0 
Ps3° ° 
P = e e s e s se 
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The weight matrix consists of weighted observations 


of a traverse with mixed kinds of measurements. They are 


distances and angles. The variance of a measured distance 


is expressed in meters and an angle in radians. For 


uncorrelated observations of unequal precision, the funda- 


mental condition of least squares is expressei as 


M 

3 
< 

—"n 
i" 


P,, 7; + Poo V5 + piv? ¢ . . « « + Pp 


Mininun 


or in matrix forn 


v'py = Mininun (3.14) 


In the problems involving the adjustment of observed 


values, all observed quantities are expressed by functions 


of the quantities to be determined. In simple cases these 


relations are linear, but when this is not the case, the 


relations must be converted into the linear form by 


expanding them into Taylor's series. The terms of higher 


order are neglected, so as to obtain linear relations, 


solving the resulting linear equations, then iterating until 


the effect of the neglected higher order terms are 


Minimnized. 


2. Least 


ice 
wm 
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The observed values are related to the jesired 


unknown values through formulas or functions which are 


called observation equations. One observation equation is 


written for each measurement. To solve for the best value 


of each unknown parameter, at least one redundant observa- 
tion equation must be written. That is, the number observa- 
tions must be greater than the number of unknowns. The 
linear okservation equations can be written in the general 
form as follows: 


t 


’ , ae 
. re ese a 6 
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a G, + ¥, - 
, G5 + Vo L) 
e e 
e e d 
Xm = G, i Vn 


where n is the number of observations; 


gm is the number of 


the unknowns; A190" Aoqe + © © Agog are constants; a,,, 4,5, 


+ + + Anm are coefficients of the unknowns X,, X5, ++ - 


X 3; and Vi « Yor + + + Vy are the residuals. 


m 


Because the observations G; (i = 1, 2, . - - n) are ; 
not free from random errors, each G; must be corrected by a 
residual value, v;- let be = G - aioe Thus 
Bouma Pasha: © oes Soin Sey = har ee 
On1%, * Oe or ie ee eda DY | 
or in the matrix forn 
voe= AX - B (3.15) 
This equation is called the observation equation or : ] 
the observation equation matrix. 
For uncorrelated observations of unegual precision, 
substitute the value for the V matrix from the observation 
Equation (3.15) into Equation (3.14) bd. 
32 
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v'pv = (AX - B)' P(AX - B) 
= [(ax)’ - Bl yp(ax - B) 
= (xTa’ - B’yp(ax - B) 
= (x’a’p - Blpy (ax - B) 
= x’al pax - x’a’pp - Bl pax + B’ PB 
= x'al pax - x’a’ pp - xTa’pp + BTPB 
= x’al pax - 2x’aTpB + BlpB (3. 16) 


[from matrix algebra, (Ax)’ = x’a? ana Blpax = x'a™ pp] 

The minimum of this function can be found by taking 
the partial derivatives of the function with respect to each 
unknown variable (i.e., with respect to the Xie Xoe eo 2 
Xm) must equal zero. Hence, 


ven) = 2a’ pax - 2a'PB = 0 (3.17) 
x 


Dividing Equation (3.17) by 2, the following result 
is obtained: 


t 


a’pax = a pp 


This represents the normal equations, and by multi- 


plying this equation by (ATPA)-1, the solution is obtained: 
xX = (aA’pay—1taTpB (3. 18) 


For uncorrelated observations with equal precision, 
the weight matrix is an identity matrix, and equation (3.18) 
becomes 
x = (al ay-1alB (3. 19) 


This equation can be derived similarly to the 
unequal precision case. Fqguations (3.18) and (3.19) are the 
tasic least squares matrix equations. 
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If the relations are nonlinear, the relations must 
be converted into the linear form by using Taylor's series 
expansion [Ref. 3, p. 919]. Let F=f (X,- Xone ee i) 
be the general observation equation that is a nonlinear 
function. The Taylor's series expansion is 


Fs £ (x0 X50 6 eX) + Of bx,+ OF See a oe OE O x 


2° 
Ox4 0x, Ox", 


m 


+ Higher order terms 


where X\, X5, - + - Xm are the approximate values of the 
variables at which the function is evaluated. For the 
traverse problems, an approximate value can be the 
precalculated value, and &x,, One Ss verve Sx, can be the 
corrections. The higher order terms in the series are 
neglected and only the zero and first order terms are 
Maintained, the approximate value must be improved by 
successive iterations until the erfect of the neglected 
higher order terms is minimized. After linearization, the 
observation equations become: 


WS 2p (kos Zoe eos 58) +9£,6x,+ 3£,5x,+ - - +O£,5x,- G, 
° x? 


xy i) Ye) m 
v, = £,(x)6 x0 2-2 Xt) + Of,0x,+ 2£,5x,+ ie ose +O¢,5x - G, 
Ox, Ox‘, Ox 


Vv, = f(t) 0 Xe 2 XH) +£,5x,+ Of,5x,+ - . +Of Sx -G 


n Pas ° ° n m nt 
Ox’, Ox’, Ox, 
or in the matrix forn 
vVo= AX- 8B (3.20) 
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The remainder of the least squares procedure is the 
same as indicated by Equations (3.16) to (3.19). 

There are two types of observation equations for the 
adjustment of the traverse. They are the angle and distance 
observation equations. Both of them are nonlinear, and they 
must be linearized using the Taylor's series. 

To derive the observation equation for the angle a, 
(Figure 3.1), first write the two azinuth observation 
equations for the azimuth Az, and AZi: , where AZ is the 
back azimuth and AZ; is the forward azimuth. 


Ea Ey 
= pas ieee 
AZ: arc tan = N; 
Ej Ses ceca 
Az.. = arc tan Ny NN sae 
i : 


ij 


wai te ee a a Pl 
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Then, the angle observation equation for the angle a; is 
v. = AZ;; - AZ a. mae 


» . —E-; 
arc tan eo - arc tan &kK-~i  - at, (3.21) 
NIN; N,—N; i 


Fgquation (3.21) is nonlinear in the parameters. 


Write this equation in function forn: 


vy, = £; (Bye Nye Eye Nin Eye NY) - 
Thus, the linearized form of the angle observation 
equation is ‘ 
vy, = £: (Ee Nps Bp, Ne Ee Ny.) +E bz, +2£,5N 
i ONy 
+Of bE, + Of, Oy, +2¢,5e, +3f,5n, - a: oe 
OE ON? E° N;, . 
The observation equation for the distance d; between ie oe 
two points i and j is < 
ra 
v,; = ((E, - E;)? + (N; = hae aed =a - : 
Fae 
or ir. the functional form ay, 
Vi FE (Epo Ngo Eye Ny ) - 4; oe 
Thus, the linearized form of the distance 
observation eguation is 
v, = £;(B., Nie Bi, NW, ) + O£-58, + O£,6n, 
J J 3 | ! > 3! | 
+ Of,58;+ Oty, - a, 
OE* Ons 
J J 
® 


For the adjustment of the traverse which was 
conducted at Moss Landing with the technique of adjustment 


of indirect observations, four angle and three distance 


< observation eguations need to be written. The seven equa- 


tions would include as unknown parameters the four coordi- 
nates of stations 2 (Mossback) and 3 (Dune Temp). Four 
normal equations are formed and solved for corrections to 
the approximate values (precalculated values) of the parane- 
tecs. The corrections are added to the approximations to 
update their values, and the solution is repeated until the 
last set of corrections is insignificantly small. 


The linearized eguations of this problems are 


v= a,,5E, + a,,5N, + a,3 SE, + ay, by, - b, 
5 a,, OE, + a,, oN, + a,3 SB, + a, , ON, - b, 
v= a,,0E, . a,,0N, * a3 bE, +a,,5N, - by 


“The zero order terms of the angle and distance 


functions are 


E,-E 
= egies el, et 
f, arc tan N5—Ny AZ. 
[ree oe <I ° 
f, = arc tan —#——2 - arc tan £17 £2 
2 N3- No Ny—No 
3 = Ee E3— E* 
f. = arc tan ~4——2 - arc tan —2___3. 
3 N N NS — N 
ao ae ae 
E3-E, 
f. = AZ, - are tan yey 
3 4 
f£, = ((E, - EB }2 (8, = hs ets 
f, = ( (E, _ EB’, )2 + (Ne, eae ee baie 


rh 
' 


7 = CE, - B, 2 + (Wy - Np 2772 


By using the data of the known coordinates (Table II) and 


the approximate values (Table VII) 


a,, =f, = N> - Ny = 0.68246 x 10-3 
OE ONS — NW)? + (ESE, )* 
2 1 2 1 
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® 
aj, =f, = - E> - £; i= 9.16550 x 10-3 
ONS (NZ-N,)? + (E5- £,)? . 
a43 “28 = 0 
3 
a4 = 28 teu ® 
3 
ao, =2f, = - N3 - NS + ae ae 
Oe (N3-N3)° + (E4- ES - (Nj- NG 2 + (E,— £5)? 
) 
= - 3.06868 x 10-3 
ee ae E3 - E> a E, - E> 
N° _ 2 ae) _ ae \2 _ ee 2 
2 (No N5 ) + (Es E,) (N, N,) + (E, E 4) = 
= 1.16926 x 10-3 
a4 S063" = uy - ND = 2.38621 x 10-3 
23 2 i+... “4 
GES KN3- NB)? + (E5- E53)? ne 
» 
an, =2f, = - E3 - E> = - 1.33476 x 10-3 
Ong (N2-N3)? + (E° es 
a3, =0f, = - Nz = 2.38621 x 10-3 Bo 
ES (N53 -N3)? + (E8 ae. 4 
azo =2f; = E2 - £3 = - 1.33476 x 10-3 
ONy (No-N%? + (E5- £3)? 
_ _ ° Z $ : ) 
433 = f5 = “f Na = N3 + N>5 - N3 
E3 (N= NS 2 (EB - £42 [NB - N3)? + 5-85) 
= - 2.50541 x 10-3 es 
=7% = ° a) ° r) , 
a4 “3B = | Eq - E3 i] | E> - £3 
2. ine 2 ° °o 4f \2 ° ° \2 
3 (Ny—NQ)® + (E,- £594) NNR He) 
= 1.43291 x 10-3 
® 
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aes = eed ° a = Pe 

a, = oF, sun Me 5 0.11920 x 10 
| 3 (N3- N,)%+ (ES-E4) 


a,, =f Er B3 - Ea ic - 0.09816 x 10-3 
( 


a,, =f, = | E1_- 22 i - 235.67466 x 10-3 
° \2 TNRAD ] Ive 

ED (e ,- €)? + (N, —N%)? ] 
a., =3f5 = - N) - No = 971.83201 x 10-3 
Ons o 2 ° 2\72 

(E,- E5)° +(N, — NS) E 

ag, =2f5 = 0 

EA 
a., = of. = 0 aoe 
54 One 


a2) Ol. Ep - £3 488.17948 x 10-3 


° ° (2 ° ° 2 
lles-e%,) +(No-N3 ) 


| a ) 
= 
‘\ 
wily 
" 
' 


872.74326 x 10-3 


a., =Ofg = - E> - E° = 488.17948 x 10-3 
ne rol = ° 2 ° 2|172 
3 (E5-E3)° +(NZ—N4) 
a_, = 9fg = - N> - N53 = 872.74326 x 10-3 
of ONS see e 2 ine he y2 17/2 
Seg Eg rea NS 


a =of,=0 


o 
D 
— 
' 
lo) Foy) 
4 lm 
0 
N 
= ’ 
m 
ho 
i 
mM 
Wo 
are) 
+ }>° 
za 
Nh 2 
; |= 
z 
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RS 
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E3 - Ey }= - 635.68254 x 10-3 
lene? +as_—n, Pe ] 4 
Ny - N - 771.95059 x 10-3 
One (eo o ee 
3 fex-e,% + 0NS =n) 


or in the matrix forn 


0.68246 0.16550 0.00000 0. 00000 
-3. 06 868 1.16926 2.38621 -1.33476 
2.38621 -1.33476 -2.50541 1.43291 
0.00000 0.00000 0.11920 -0.098 16 
-235.67466 971.83201 0.00000 0.00000 
-448.17948 -872.74326 488.17948 872.743 26 
0.00000 0.00000 -635.68254 -771.95059 


The computation of B matrix, by using the data 
from Tables II, III, IV, and VII 


b, 0.00000 radians 
0.00000 radians 

0.00002 radians 

0.00002 radians 

0.00000 meters 

0.00000 meters 

-0.01426 meters 


or ir the matrix form 


B' =[{ 0 0 0.00002 0.00002 0 0 -0.01426 } 


The weight matrix can be computed by using the vari- 


ance of the observed angles and distances from Tables III 
and IV by applying to Equation (3.13). For uncorrelated 
observations of unequal precision, the weight matrix is 
defined as the diagonal P matrix. 
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PR, = 17 sin2(1.984") = 10,808,537, 426.79957 
Poo = 17 sin2(1.405") = 21,552,498,219.02474 : 
Pj, = 1 7 sin?(1.203") = 29,398,082,997.43487 | 
P,, = 17 sin2(1.619") = 16,332,144, 194. 08730 
Pa,= 1 7 (0.001)2 = 1,000,000.00000 
P= 17 (0.0012 = 1,000,000. 00000 ; 
p..= 1/7 (0.003)2 = 111, 117.11111 


For uncorrelated observations of unegual precision, 
the values of A, B, and P matrices are applied to Equation 
(3.18). The correction vector is found to be 


SE, = 0.01284 
xX = 5x, = 0.00336 
Sz, = 0.01079 
SN, = 0.00495 


The approximation values were improved by adding the 
corrections to the first approximation values. The improved w 
approximation values after the first iteration are oe 


E, = 607,943.44238 + 0.01284 
5 4,073,939.74473 + 0.00336 
+ 
+ 


607,943.45522 m. 
4,073,939.74809 m. 
608,122.00189 m. 
4,074,258.95029 n. 


=z 
It 


608, 121.99110 0.01079 
3 4,074,258.94534 0.00495 


= bf 
WwW 
i 


Using these values, the solution is iterated. After the 
second iteration, the correction vector is zero to six 
decimal places, so the improved approximation values are 
the final estimates of the coordinates. 

For uncorrelated observations of equal precision, 
the correction vector is found by solving Equation (3.19). 
The correction vector after the first iteration is 


SE, = 0.01756 
X= BN, = 0.00426 
SE, = 0.01660 
Su, = 0.00479 
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The impreved approximation values after the first iteration 


IES <n, Ol , ae Ro 7 
| 


are ‘ - 
E. = 607,943.44238 + 0.01756 = 607,943.45994% am. 
N. = 4,073,939.74473 + 0.00426 = 4,073,939.74899 a. ore 
E, = 608, 121.99110 + 0.01660 = 608,122.00770 am. we 
N, = 4,074,258.94534 + 0.00479 = 4,074,258.95013 a. 


The solution is iterated ty using these adjusted 


values. The correction vector is zero to six decimal places 
after the second iteration. These values are the final 
estimates of the coordinates. 

The standard deviation of an observation which has 
unit weight can be found by the following equations, for 


uncorrelated observations of unequal precision, ‘ 


6, = [ er] 172 (3.22) s 
nA-m 2 
for uncorrelated observations of equal precision, ae 
6, = [ee] 172 (3.23) a 
n-m 
re: where n is the number of observation eyuations and m is the 


number of unknowns [Ref. 7, p. 249]. After calculating the 
best estimate values of the unknowns, the V matrix can be 
computed from Equation (3.20). 

The standard deviations of the best estimate values 
for the unknowns are then given by the following equations: 


6, = 6,[s::] 1/2 (3. 24) ot 


where 6, is the standard deviation of the ith adjusted 
quantity. The quantity in the ith row of the X matrix, S.. 
is an element of the (A’ PA) -! matrix for uncorrelated 2 
observations of unequal precision. For uncorrelated 

observations of egual precision, S is an element of the 


(ATA) 1 matrix [Ref. 7, p- 250]. 


oe 
Vt 
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The standard deviations of adjusted angles and 


distances can be computed by the following equation: 


6. = 6.[ 243] 172 (3.25) 


where 6, is the standard deviations of the ith adjusted 
guantity. The guantity in the ith row of the V matrix, 9: 
is an element of A(A’ pay -tar matrix for uncorrelated 
observations of unequal precision. For uncorrelated 
observations of equal precision, Q,, is an element of the 


A(A’ A)—1a" matrix [Ref. 3, p. 912]. 


A(aA’ PA) -2aT é J Be See Gee gh 
or = s e ee e @ oe eo e 
A(A’ A) -2aT és - bre, areal oe 


Qny Qn. . e ee ° Q 


A closed traverse was conducted at Moss Landing. 
After the second iteration, the V matrix can be found by 
applying the values of A, B, and X matrices to FEy\uation 
(3.20). 
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For uncorrelated observations of uneyual precision, 


the V matrix is found to be 


¥, = 0.93149 x 10-S radians or 1.92134" 

Vzo = 71663120 x 10-5 radians or -3.36459" cae 

V3 = —1-28363 x 10-5 radians or -2.64768" i 
voe= Y4 = 2230836 x 10-5 radians or -4.75308" 

Ya. = 0.00024 meters 

Yo = 0.00039 meters 

YW = 0.00358 meters 


The standard deviation of unit weight is found by solving 
Equation (3.22), the result is 


6, = + 2.69685 
The values of cal PA)—2 are ' 


1.25178 0.22019 1.25575 0.10788 4 
10-5 | 0.22019 0.13819 0.21810 0.11127 ae 

1.25575 0.21810 1.96408 0.03815 

0.10788 0.11127 0.03815 0.22515 


The standard deviations of the best estimate values for ' 
positions can be found by solving Equation (3.24), the 
results are 


6-5 = + 0.00954 seters 
Oyo = + 0.00317 meters =: 
E3 = + 0.01032 meters oe 
N3 = £ 0.00405 meters 


4y 


The values of A(a’pA)-1Aa" are 


0.064 -0.052 -0.021 0.010 -4.165 -6.963 -63.88 

“0.052 0.225 -0.175 9.002 3.415 4.690 41.59 

-0.9021 -0.175 0.211 -0.914 2.563 3.574 31.79 

19-10 9.010 0.002 -0.014 0.002 -1.813 -1.301 -9.501 
“4.165 3.415 2.563 -1.813 9917 -117.4 -1047 

-6.963 4.690 3.574 -1.301 -117.4 9830 +1525 

-63.88 41.59 31.79 -9.501 -1047 -1525 76323 


The standard jeviations of adjustei angles and iistances are 


founG by solving Equation (3.25), the results are 


6 = 


a4 t 1.40346 seconds 
6,5 = + 2.64103 seconds 
6,3 = + 2.55260 seconds 
O34 = + 0.26137 seconds 
O51 = + 0.00269 meters 
Sy. = t 0.00267 meters 
O53 = + 0.00745 meters 


For uncorrelated observations of equal precision, 
the V matrix is found to be 


vi = 1.26896 x 10-S radians or 2.61742" 

Yao = 71256893 x 10-% radians or -3.23615" 

Va3 = 71275914 x 10-5 radians or -3.61817" 
voe= V4 = 72.23358 x 10-5 radians or -4.60710" 

¥, * 0.149037 x 10-7 meters 

es 0.23844 x 10-7 meters 

A a 0.24365 x 10-7 meters 


The standard deviation of unit weight is found by 


solving Equation (3.23), the result is 


: O, =+ 2.01144 x 10-5 


The values of (av A) 2 are 


0.28193 0.06827 -0.85449 0.70380 

10-4 0.06827 0.01664 -0.20731 0.17075 
-0.85449 -0.20731 2.59561 -2.13724 
0.70380 0.17075 -2.13724 1.75999 


The standard deviations of the best estimate values for 
positions can be found by solving Equation (3.24), the 
results are 


O-> = + 0.00107 meters 
Gyo = + 0.00026 meters 
Oc3 = + 0.00324 meters 
Gis = + 0.00267 meters 


The values of A(a’ ay—tal are 


[ 1,472 -27.197 26.960 -1.235 -0.518 -0.343 -0.356 
-27.197 503.1 -498.7 22.855 9.348 0.480 0.473 
26.960 -498.7 494.4 -22.66 0.370 0.523 0.518 
10-3 “1.235 22.855 -22.66 1.038 -0.201 -9.159 -0.134 
-0.518 0.348 0.370 -0.201 10000 -0.001 -0.001 
-0.843 0.480 0.523 -0.159 -0.001 10000 -0.001 
-0.856 0.473 0.518 -0.134 -9.001 -0.001 10000 


The standard deviations of adjusted angles and distances are 


found by solving Equation (3.25), the results are 


O21 = + 0.15917 seconds 
O32 = + 2.94270 seconds 
©.3 = + 2.91732 seconds 
6. = + 0.13370 seconds 
orn = + 0.00002 meters 
Sue = + 0.00002 meters 
O33 = + 0.00002 meters 


3S a 4 Hoe ag aes 
Pee Os pane ye 
ee er ee ce ae! Lae at rd Pe 


The computations in this subsection were performed 
by the NPS IBS 3033 computer using 16 decimal places and 
were rounded off to 5 decimal places prior to output. 


3. Least Squares Adjustment by the Condition Equation 


————— 


The general principle of least squares adjustment by 


the condition equation method in surveying is 


to minimize a function consisting of the sum of the 
Squares of the corrections to the observations plus the 
necessary mathematical conditions involving some or all 
the corrections; each condition by itself is nade egual 
to zero by adding the corrections to the discrepancy 
determined from a preadjusted calculation (i.e. calculated 
from the observed values, rather than from adjusted ; 
values); thus, the magnitude of the sum of the sSguares is 
not changed by adding conditions {Ref. 8]. 


For uncorrelated observations of unegual precision, 
the principle condition of least syuares function may be 
expressed in matrix form as follows: 


U = Vi pv - 2K(BV + W] 


= Vipy 


2[BV + WI'K 


= v'pv - 2¢(Bv)’ + wT jK 


= vi py - 2cv’BT+ w' yK 
= v’pv - 2v'’B’K + 27K 
= mininun 


where U is the least squares function matrix, K is a matrix 
of Lagrange multipliers, P is the weight matrix, V is the 
correction vector, B is the constant coefficients of the 
corrections, and W is the precalculated discrepancy 

(Ref. 9]. It is numerically more convenient for later 
development to multiply by 2. Taking the partial 
derivatives of the 0 matrix with respect to each of the 
corrections and eguating to zero leads to 
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Qu = 2Pv - 2B)K = 0 
Ov 
or 
Vv = Pp-1B'K (3.26) 


This equation is called a correlate ecuation or correlate 
eguation matrix. The solution of the Lagrange multipliers 
matrix can be derived by multiplying Equation (3.26) by the 
B matrix then adding the ¥W matrix 


BV ¢ W = BP-iB'K + W 


(BV + W) is the condition equation matrix which must egual 


zero. The solution of the Lagrange multiplier vector is 
BP-1B’ K = -W 


K 


(BP-1BT) -1 (-8) (3. 27) 


The correction vector is derived by substituting 
Equation (3.27) into Equation (3.26) 


Vv = p-1B7 (BP-1BT)-2 (-§) (3.28) 


For uncorrelated observations of equal precision, 
the correction vector can be derived similarly to the 
unequal precision case or by using the inverse of the weight 
Matrix which is the identity matrix. Equation (3.28) 


becomes 
v = B! (BBT)-1 (-w) (3.29) 


Closed traverse station positions may be adjusted by 
using the technique of least squares adjustment by the 
Condition Equation Method. There are two condition equa- 
tions. They are the azimuth and coordinate condition equa- 
tions. The coordinate condition ejuations are divided into 
two parts, which are E and N coordinate condition equations. 
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The azimuth condition equation is the sum of the 
corrections to the angles in a traverse plus a precalculated 
discrepancy W, (Equation 3.6) which must equal zero. Where 
Vai 
tion can be written: 


3. en ee 


21 
n 
2h + 20-44; + ape (3-390) 
ts 


equals the corrections to the angles, the general equa- 


or ; 


where Vgj is the correction to the distances. However, 
this term equals zero and does not effect the equation, but 
it does reserve space in matrix notation for distance 

} corrections. A closed traverse which was conducted at “oss 


Landing can be written with an aziauth condition of 


1.V + lev + ev + T.v + O.V + 0.V + O.v, 


d2 3 


] + W, = 0 (3-31) 
where the precalculated discrepancy (W,) is +8.844" or 


+0.00004288 radians. 
The simpic linearized form of an E coordinate condi- 


tion equation is the sum of the adjusted departures and rust 
equal zero. By applying Equations (3.2) and (3.7), the 
general formula for E coordinate condition equation is 
n-1 
) 2(4; + ¥,;) sin Az, - (2, - E,) = 0 (3.32) 
where —~ 
Az; = Az, + DA, + v,;) - (i - 1) 180° (3.33) 


jz 
’ Az,; is the adjusted azimuth. By using Equations (3.32) and 


(3.33) with the data conducted at Moss Landing, these 
equations become 
(dq; + ¥y,)-sin Az, + (d, + ¥,5)-sin Az,, 


- E,) = 0 (3.34) 


> + (d, + ¥43)-Ssin AZ.3 7 (EL 
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where 


Az., = Az, + A, + v4 
= AZ, + Va; 
AZ, = AZ + a, ae a, + V,. ~ 180° 
i AZ, id Vat - Va2 
AZ,, = Aze + O, + ¥,, + A, + v,, + A,+ v,, - 3600 


Az, is the precalculated azimuth or unadjusted 


azimuth (Equation 3.1). 
From trigonometry 


Sin (A+tQOA) = sin A.cosQA + cos A.sinfAA 
in which AA is a very small angle in radians, let 
cosMA = 1 and sinQaA =Oa 


therefore, 
sin (A+MOA) = sin A +QA.cos A 
from the first term of Equation (3.34) 


(d, + Vy,)Sin Az,, 


(d, + Va Ysin (Az, + Vay) 


= (d, + Vg, (Sin Az, + V4, C08 Az, ] 


d,sin Az, + v,, d,cos Az, + v,, 


sin Az, + Wj, ¥,, Cos Az, 
(3.35) 
(Vo, -V51) is very small, by letting this term equal zero, 


Equation (3.35) becomes 
(d, + Vy,) Sin Az, 


= d,sin Az, + V,,4,cos Az, + Vj, Sin Az, (3.36) 
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The second and third terms of Equation (3.34) can be derived 


| <i I 


similar to Equation (3.36), so ’ 
(4, + Vy>)Ssin AZ,> 

= d,sin Az, * (v,, + V2) 4, cos AZ, + V,, Sin AZ, (3.37) 

(d., + V,,)sin Az., 


= d,sin Az, + (v,, *¥,> +V, 3) d,cos AZ, + V3 sin Az, (3.38) 


L Substituting Equations (3.36), (3.37), and (3.38) 


into Equation (3.34) and rearranging terms 


Vai (d,cos Az, + d,cos Az, + d,cos Az.) 


- + vi (d,cos Az, + d,cos Az) + V,34300S Az, = 
+ Vy, 548 Az, + Vqo San AZ, + V43 Sin AZ, 
iH + d,sin Az, + d,sin AZ, + d,sin AZ, - (E, - E,) = 0 
- 
4 
This equation can be written in the gereral formula as 
n-1 n-1 -1 n-1 
> (V5; d,cos Az,) # Dy; sin Az; + >4,.sin Az. 
a iz kz i 121 i:1 
- (8, - £,) = 0 (3.39) os 
substituting Equations (3.2), (3.3), and (3.7) into Fquation 
(3.39) 
, n-1 n-1 n-1 2 
>(¥,; >ON,) + Dv, sin Az, + W, = 0 (3.40) 
| ket j=1 
> For example, if the number of the observed angles (n) at ; 
: Moss Landing is four, this equation can be expanded to = 
Vaq (ON, + OND + OR) ce te Gas +ON,) + v,,0N, 
: + V4,S51n Az, + Vy Sin Az, + Vy3 Sia Az, + Wo = 0 
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Substituting the numerical values from the precalculated of 
this traverse into this equation, the result is 


(6702.37612)¥,, +(5318.48345) v,5 +(4999.28284)¥ , + 0.v., 


-(0.23567)v,, #(0.48818)¥,, +(0-63570)v¥,, # ¥, = 0 (3.41) 


where the precalculated discrepancy (W,) is 0.08608 m and 
O.v,,is equal to zero, which does not effect this equation, 
but space is reserved in matrix notation for the last angle 


correction. 


The N coordinate condition equation is the sum of 
the adjusted latitudes and must equal zero. By applying 
Eguations (3.3) and (3.7), the general condition equation is 


n-1 
> (d; + vj) cos Az,, - (N, - Ny) = 0 
i:1 


This equation can be derived similarly to the £ coordinate 
aa oer pao: resulting in 


Simi 208) + o ;cos Az, + W, = 0 (3.42), 


i: kz1 
From the data Peano at Moss Lanijiing, this eguation can 


ke expanded to 

-¥,1 (OE, +OE, + OE) - vz2 (OE, +OE,) + v,3 OF, 

+ ¥41 COS Az, + Yyo cos AZ, + Vaco= AZ, + W3, = 0 
Substituting the numerical values from the calculation of 


this traverse into this equation, the result is 


~ (3959.89460)v,, -(4295.49626) v,, -(4116.94755)v., + O.v., 


+(0-97183)v,, + (0.87274) ¥,, +(0.77194) v4, # HW, = 0 (3.43) 


where the precalculated discrepancy (W4) is -0.08936 nm. 
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Equations (3.30), 


coefficients of the corrections, 


(3.40), and (3.42) can be written 
in the matrix form as (BV + W), where B is the constant 


Vis the correction 


vector 


which is composed of the angle and distance corrections, fW 


and (3.43) are written in 


6702.37612 
5318. 48345 
4999.28284 
0.00000 
-0.23567 
0.48818 
0.63570 


[20 CO = = wt mt 


' 
2 
iT] 


; = ~0.00004 
= ~0.08608 
0. 08936 


i] 
q 
It 
‘ ' 
qQ x 
"ot 


is the precalculated discrepancy. Equations 


Matrix form as 


-3959.89460 
-~4295.49626 
-4116.94755 
0.00900 
0.97183 
0.87274 
0.77194 


(3.31), 


(3.41), 


For uncorrelated observations of unequal precision, 


the inverse of the weight matrix is the diagonal matrix, 


which is equal to the observation's variance. 


The diagonal 


elements of the inverse of the weight matrix is 


x 


x 


x 


10-10 


10-10 


10-10 


10-5 


P7 = sin2(1.984") = 0.92519 x 10-10 
P = sin2(1.405") = 0.46398 
bt = sin2(1.203") = 0.34016 
Pe = sin2(1.674") = 0.61229 
-1 = (0.001) 2 = 0.10000 
Pe ( ) 
p-12 = (0.001)2 = 0.10000 
66 
-1 = (0.003)2 = 0.90000 
Pes ( ) 


ee a eg eet at Ue a er ea aad et a 
ee ae le oy Om BO OT ee wee ee ig Oe ate, 28 
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For uncorrelated observations of unequal precision, 


the correction vector is found by solving Equation (3.28) 


Va, = 0.93741 x 10-5 radians or 1.92117" 

Vao = 71263114 x 10-5 radians or -3.36446" 

Va3 = 71.28358 x 10-5 radians or -2.64758" ie 
ve Ya4 = 72.30438 x 10-5 radians or -4.75312" a 

¥,, = 0.00024 meters 

Yuo = 0.00039 meters 

Vg3 = 0.00357 meters 


=) 
tr 
© 


corrected angles are ’ 


QO, = 2460 05" 45.121178 
QO, = 2220 51° 05.23554" 
A, = 1900 14" 59.95242" wis, 
Ch, = 2779 05" 12.24688" 


The corrected distances are 


1424.00424 am. ae 


d, = : 
dp = 365.74439 ao. 
d; = 6476.27457 2. 


Using the corrected data to recompute the coordinates is 


Similar to the initial traverse computations. The adjusted 


coordinates are 


E> = 607,943.45521 a. 
KR, = 4,073,939.74809 om. 
E, = 608,122.00189 am. ' 
HR, = 4,074,258.95029 am. 
E, = 612,238.85256 am. 
BH, = 4,079,258.31754 om. 


The last adjusted coordinates are the same values as the 
fixed coordinates. Hence, it is not necessary to iterate. 
For uncorrelated observations with equal precision, the 
correction vector is found by solving Eguation (3.29) 
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The corrected angles 
a, 
Xs 


a, 
A, 


1.26884 
-1.56886 
-1.75408 
-2. 23360 

0.14034 

0.23842 

0.24362 


~* « Mw HM OO 


= 2460 05! 
= 2229 51! 
= 1900 148 
= 2779 05° 


The corrected distances are 


d 


a 


By using the corrected 


the adjusted coordinates are 


Eo 
Ny 
Ey 
M3 
E4 
Ny 


The last adjusted coordinates are the same values as the 


fixed coordinates. 


Hence, it is not necessary to iterate. 


data to recompute 


10-5 radians 
10-5 radians 
10-5 radians 
10-5 radians 
10-7 meters 
10-7 meters 
10-7 meters 


45.81717" 
05. 36399" 
58.98196" 
12.39287" 


1424.00400 ao. 
365. 74400° m. 
6476.27100 am. 


607,943.45994 ao. 
4,073,939.74899 nm. 
608,122.00770 am. 
4,074,258.95013 ao. 
612,238.85256 ao. 
4,079,258.31754 a. 


oc 
or 
or 
or 


the coorijiinates, 


2.617178 --- 
- 3.23601" 

-3.61804" 
-4.60713" 


The standard deviation of an observation which has hi 
unit weight can be found by the following equations. For eae 
uncorrelated observations of unegual precision,  ] 


6, = [eee] +72 (3.44) 
r 
and, for uncorrelated observations of egual precision, , 
6, = [ata] 72 (3.45) 
r i 
where r is the number of condition eguations. There are . 


three condition equations for a closed traverse (r=3). The 
standard deviations of adjusted angles and distances can be 
found by the following equation: 


6. = 6.{ 0,4] 172 (3.46) 


where 6; is the standard deviation of the ith adjusted = 
guantity. The quantity in the ith row of the V matrix, Q 
is an element of the P-t - p-1B'(pp-1B")-1pp-1 matrix for 


uncorrelated observations of unequal precision 


- 
uy 


(Ref. 3, p.- 918]. For uncorrelated observations of equal 
precision, Qj; is the element of the I - B’ (BBT) -18 matrix. . 
The matrix P-1 - p-1B) (Bp—1B1)—13p-1 or the matrix 
i.e al (BET)—1B is equal to 


S25 See are 


QO Qt2 2 ee ee Qn 
Car Ogg tse Soe Ron 


e e e e e e e e 
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For uncorrelated observations of unegual precision, 
the standard deviation of unit weight is round by solving 
Equation (3.44), the result is 


6, = + 2.69679 


The results of P-! - p-13) (pp-137)-1Bp-1 are 


0.064 -0.052 -0.021 0.910 -4.165 -6.963 -63.88 

-0.052 0.225 -0.175 0.902 3.415 4.689 41.59 

-0.021 -0.175 0.21171 -0.014 2.563 3.574 31.79 

19-10 0.010 0.002 -0.074 0.002 -1.813 -1.301 -9.500 
-4.165 3.415 2.563 -1.813 9917 -117.4 -1047 

-6.963 4.689 3.574 -1.30171 -117.4 9830 -1525 

-63.88 41.59 31.79 -9.500 -1047 -1525 76323 


The standard deviations of adjusted angles and distances are 
found by solving Equation (3.46); the results are 


al +t 1.40340 seconds 
6.5 = + 2.64096 seconds 
Os3 = + 2.55253 seconds 
6,4 = + 0.26136 seconds 
Og, = + 0.00269 meters 
Oa = + 0.00267 meters 
643 = + 0.007485 meters 


For uncorrelated observations of equal precision, the 
standard deviation of unit weight is found by solving 
Equation (3.45), the result is 


O, = + 2.01139 x 10-5 
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The results of I - BY (BBY y—1B are 


J 
1.472 -27.199 26.962 -1.235 -0.518 -0.843 -0.856 a 
-27.199 503.1 -498.7 22.855 0.348 0.480 0.473 oe 
26.962 -498.7 494.4 -22.66 0.370 0.523 0.518] re 
10-3 -1.235 22.855 -22.66 1.038 -0.201 -0.159 -0.134 
-0.518 0.348 0.370 -0.201 10000 -0.001 -0.091 
-0.843 0.480 0.523 -0.159 -0.001 10000 -0.001 
-0.856 0.473 0.518 -0.134 -0.001 -0.001 10009 ; 
The standard deviations of adjusted angles and distances are 
found by solving Equation (3.46); the results are ey 
: 
O,; = + 0.15918 seconds ee 
O,> = + 2.94262 seconds 
6,3 = + 2.91722 seconds _ Tes 
(G4 = + 0.13369 seconds — 
O,, = + 0.00002 meters . 
66 = + 0.00002 meters ite 
643 = + 0.00002 meters Te 
The computations in this subsection were also > 


performed on the NPS IBM 3033 computer; they used 16 decimal 
places, and were rounded off to 5 decimal places in the 
final solutions. These are more decimal places than are 
normally used in practice because such precisions are not 4 
generally attainable ty corresponding observations. 

However, in comparing two computational methods that are 

theoretically equal, one method may be more sensitive to 

rouné off error than the other. This will be commented on . ® 
in the next chapter. 
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Three programs were used to compute the initial traverse 
and adjust the closed traverse station positions. Programs 
1, 2, and 3 were used to compute and ajijust the traverse 
station positions by the Approximate Method, the Indirect 
Observations Method, and the Condition Equation Method, 
respectively. The programs were written in WATFIV language 
for implementation on the NPS IBM 3033 computer. The 
computer output and listings of Programs 1, 2, and 3 are 
provided in Appendices A, B, and C, respectively. The 
waximum number of intermediate traverse station positions 
that can be computed and adjusted by these programs is 30. 
The programs were tested using several fictitious data sets 
to ensure their performance in handling the various 
intermediate traverse station positions. 

The computer storage area and CPU time of Programs 1, 2, 


and 3 has been compared (Tabie X). The adjustment of a 


Closed traverse by the Approximate Method did not use a 
weight matrix in the matrix computations. The computer 
program for this method was written by using the variables 
in only one dimension for storage of both data and results. 
This method did not require ary iterations. The least 
Squares adjustment of a closed traverse characteristically 
is used to simultaneously eliminate closing errors in 
azimuths and distances, and Programs 2 and 3 -oth utilized 
such an adjustment. The computation by either Least Squares 
Method used matrices’ in two dimensions to compute the 
correction vector. Likewise the computer programs were 
written using variables in two dimensions. Included are all 
subroutines from Program 1 and extra subroutines for each 
individual program. Therefore, both Programs 2 and 3 used 
more computer storage area and CPU time than Program 1. 
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The Indirect Observations Method requires that the 
Rumber of observation equations be equal to the sum of 
observation equations of both observed angles and distances. 


i: ine nehiaiadn, vie = 
. 
294 


The number of observation equations is not constant. [It 
changes depending on the number of the traverse stations-- 
this means that the row dimension of each matrix will change ae 
; depending on the number of observation equations. ’ 
i Therefore, the subroutines affected by the number of trav- 
erse Stations are difficult to write in the general program. 
The coefficients in the A matrix are computed by taking the 
pactial derivatives of the functions with respect to each ' 
unknown variable, where the number of unknown variatles is 
not constant. The number of unknown variables changes 
dependiny on the number of the traverse stations. 
Therefore, the subroutine for the computed A matrix is also ) 
difficult to write in the general program. The Condition 
Equations Method has only three equations which makes it 
easier to derive the general form and write the computer : 
program. The adjustment using the Indirect Observations a 
Method does not apply the corrections to the observed values 
directly; therefore, it requires the approximate values of 
the unknown coordinates be computed for the correction 
vector of the unknown coordinates. This method requires at 
least two iterations to check the insignificance of the 
correction vector when it is compared with an arbitrarily 
selected small number. The adjusted coordinates from the 
Condition Equations Method can be checked at the first iter- 
ation. Thus, Program 3 is more economical as it uses less 
computer storage area and CPU time than Program 2 (Table X). 

The closed traverse at Moss Landing originated and 
terminated at control points with known positions which were 
determined by third-order methods. The azimuth closure was : 
2.211" per station and the position closure was 1:66,617. eae 
The classification, standards of accuracy, andi general wane 
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specifications for a third-order class I traverse (Table I) 
indicate the azimuth closure is not to exceed 3" per station 
and position closure must be Letter than 1:10,000. This 
traverse met the specifications and standards of accuracy 
for a third-order class I traverse. 

Using Programs 1, 2, and 3, the traverse was computed 
and adjusted in UTM grid coordinates. The final adjusted 
coordinates were transformed from UTM grid coordinates to 
geographic coordinates [Ref. 4, pp. 319-321]. The differ- 
ence between the coordinates for each method and NOS 
(Ref. 10] have been computed (Table XI). The technique of 
Least Squares for both methods yields identical computa- 
tional results to at least five decimal places (Table XI, 
Methods 3.1 and 4.1; 3.2 and 4.2). Least Squares provides 
the best estimates for positions of all traverse stations. 
Program 2 performs a statistical test, yielding the standard 
deviations of both adjusted positions and observed values. 
Conversely, Program 3 only yields the standard deviation of 
adjusted observed vaiues. Comparisons of the standard devi- 
ations of adjusted observed angles and distances were made 
between these computed by the Indirect Observations Method 
and the Condition Equations Method (Table XII). The 
computed standard deviations of the adjusted angles differed 
Significantly in the fourth decimal place in several cases, 
which is a slight indication that the Indirect Dbservations 
Method is more sensitive to round off error. The standard 
deviation estimates are larger for the Indirect Observations 
Method in every case of a significant difference in the 
fourth decimal place. However, a significant difference in 
the fourth decimal place is insignificant in terms of the 
observational precision. The standard deviations of the 
adjusted distances for both methods yield identical results 
(Table XII). Both Methods differ from the NOS positions in 
the fourth decimal place in seconds of arc for both latitude 
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and longitude. The horizontal control for third-order stan- 
dards recuires accuracy to three decimal places in seconds 
4 of arc for both latitude and longitude--thus, these methods 


can be used for computing third~order positions. 
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Ve. CONCLUSIONS AND RECOMMENDATION e 


A. CONCLUSIONS 


The computer programs developed by the author ar@ a 
contained ir this thesis have a variety of useful and prac- 
tical hvdroyraphic applications. In hydrographyv, geodetic 
rielé work glays a key role. dHorizontal control is neces- 
Sary in determining positions and, even in areas which 
appear to initially have adeuwuate control, additional 
control stations often need to be established after arrivul 
in tre Zield. These field-generated sontrol ;2.ints mist ce 
adjusted to fit within the existinjy survey net and to wini- : 
Mize ereurs in the Jieleo aeasuirceents. The errors can be 
minimized by these computer programs. Direct application of 
these programs would be enormously beneficial, since no Eo 
equivalent software exist at NPS to perform such adjustments a 
at the present time. 

In the Approximate Method, a traverse is adjusted by 
computations using a hand calculator. This method is suit- 
able for field computation and the adjusted coordinates meet ee 
the specifications and standards of accuracy for a third- 
order class I traverse. Although the accuracy from this 


method is less than the Least Squares Yetnod, it does not 
require a computer and it furnishes coordinates which can be 
checked in the field by the field party. 

The Least Squares Method provides the best estimates for 
Fositions of all traverse stations. Least Squares does 
require more of a computational effort than the Approximate 
Method. However, the accuracy of Least Sguares is better 
than the Approximate Method, and the Least Squares Method 
should be used for the final office computations. Both 
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Least Squares Methods yield identicai results, but the 

Condition Equations Method is more economical and easier to 

derive than the Indirect Observations Method. Therefore, \ 
the Condition Equations Method should be applied at NPS. 


B. RECOMMENDATION 


Studies should be continued at NPS to compare the 
econony of adjustment methods used in this thesis with other 


methods of least squares adjustment. 
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